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In the framework of ADM formalism, it is possible to find out eigenvalues of the WDW equa- 
tion with the meaning of vacuum states, i.e. cosmological constants, for f{R) theories of gravity, 
where f{R) is a generic analytic function of the Ricci curvature scalar R. The explicit calcula- 
tion is performed for a Schwarzschild metric where one-loop energy is derived by the zeta function 
regularization method and a renormalized running Aq constant is obtained. 
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, General Relativity (GR), together with Quantum Field Theory, is the major scientific achievement of last century. 
O ' It is a theory of spacetime, gravity and matter unifying these concepts in a comprehensive scheme which gives rise to 
a new conception of the Universe. However, in the last thirty years, several shortcomings came out in the Einstein 
> scheme and people began to investigate if GR is the only theory able to explain the gravitational interaction. Such 
T-H ' issues essentially spring up in Cosmology and Quantum Field Theory. In the first case, the presence of Big Bang 
singularity, flatness and horizon problems [l] led to the result that Standard Cosmological Model 0, is inadequate 
to describe the Universe at extreme regimes. On the other hand, GR is a classical theory which does not work as a 
fundamental theory, when one wants to achieve a full quantum description of spacetime (and then of gravity) . Due to 
this facts and, first of all, to the lack of a definitive Quantum Gravity theory, alternative theories of gravity have been 
, pursued in order to attempt, at least, a semi-classical scheme where GR and its positive results could be recovered. 
1 A fruitful approach has been that of Extended Theories of Gravity (ETC) which have become a sort of paradigm in 
C3 ' the study of gravitational interaction based on corrections and enlargements of the Einstein scheme. The paradigm 
C ■ . consists, essentially, in adding higher-order curvature invariants and non-minimally coupled scalar fields into dynamics 
' resulting from the effective action of Quantum Gravity [3, [3| ■ 

All these approaches are not the "/mZ/ Quantum Gravity^^ but are needed as working schemes toward it. In any 
^ ' case, they are going to furnish consistent and physically reliable results. Furthermore, every unification scheme 
^ , as Superstrings, Supergravity or Grand Unified Theories, takes into account effective actions where non-minimal 
couplings to the geometry or higher-order terms in the curvature invariants come out. Such contributions are due to 
^ ' one-loop or higher-loop corrections in the high-curvature regimes. Specifically, this scheme has been adopted in order 
to deal with quantization on curved spacetimes and the result has been that the interactions among quantum scalar 
fields and background geometry or the gravitational self-interactions yield corrective terms in the Hilbert-Einstein 
Lagrangian Q. Moreover, it has been realized that such corrective terms are inescapable if we want to obtain the 
effective action of Quantum Gravity on scales closed to the Planck length Q. 

Besides fundamental physics motivations, all these theories have acquired a huge interest in cosmology due to the 
fact that they "naturally" exhibit inflationary behaviors able to overcome the shortcomings of Standard Cosmological 
Model (based on GR). The related cosmological models seem very realistic and, several times, capable of matching 
with the observations [sl. [ol. [l0| . Furthermore, it is possible to show that, via conformal transformations, the higher- 
order and non-minimally coupled terms always correspond to Einstein gravity plus one or more than one minimally 
coupled scalar fields [ll|, [13, [H, • This feature results very interesting if we want to obtain multiple inflationary 
events since a former early stage could select "very" large-scale structures (clusters of galaxies today), while a latter 
stage could select "small" large-scale structures (galaxies today) [l^. The philosophy is that each inflationary era is 
connected with the dynamics of a scalar field. Furthermore, these extended schemes naturally could solve the problem 
of "graceful exit" bypassing the shortcomings of former inflationary models [lB| ■ 

Recently, ETC are going also to play an interesting role to describe the today observed Universe. In fact, the 
amount of good quality data of last decade has made it possible to shed new light on the effective picture of the 
Universe. Type la Supernovae (SNela) [13], anisotropics in the cosmic microwave background radiation (CMBR) p^ . 
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and matter power spectrum inferred from large galaxy surveys [19j represent the strongest evidences for a radical 
revision of the Cosmological Standard Model also at recent epochs. In particular, the concordance KCDM model 
predicts that baryons contribute only for ^ 4% of the total matter - energy budget, while the exotic cold dark matter 
(CDM) represents the bulk of the matter content (~ 25%) and the cosmological constant A plays the role of the so 
called "dark energy" (~ 70%) [l^]- Although being the best fit to a wide range of data [2l[, the ACDM model is 
severely affected by strong theoretical shortcomings [l^ that have motivated the search for alternative models [23l |. 
Dark energy models mainly rely on the implicit assumption that Einstein's General Relativity is the correct theory 
of g ravity indeed. Nevertheless, its validity on the larger astrophysical and cosmological scales has never been tested 
[24| , and it is therefore conceivable that both cosmic speed up and dark matter represent signals of a breakdown in 
our understanding of gravitation law so that one should consider the possibility that the Hilbert - Einstein Lagrangian, 
linear in the Ricci scalar R, should be generalized. Following this line of thinking, the choice of a generic function f{R) 
can be derived by matching the data and by the " economic" requirement that no exotic ingredients have to be added "'^ . 
This is the underlying philosophy of what are referred to as f{R) theories of gravity, see [2^, [H, [S^, [H, [2^, [s^, HH and 
references therein. However f{R) gravity can be encompassed in the ETG being a "minimal" extension of GR where 
(analytical) functions of Ricci scalar are taken into account. Although higher order gravity theories have received much 
attention in cosmology, since they are naturally able to give rise to the accelerating expansion (both in the late and in 
the early universe, it is possible to demonstrate that f{R) theories can also play a major role at astrophysical scales. 
In fact, modifying the gravity Lagrangian can affect the gravitational potential in the low energy limit. Provided that 
the modified potential reduces to the Newtonian one on the Solar System scale, this implication could represent an 
intriguing opportunity rather than a shortcoming for f{R) theories. In fact, a corrected gravitational potential could 
fit galaxy rotation curves without the need of dark matter (3^ . [33l . Is^ l . In addition, it is possible to work out a formal 
analogy between the corrections to the Newtonian potential and the usually adopted dark matter models. In general, 
any relativistic theory of gravitation can yield corrections to the Newton potential (see for example [36]) which, in 
the post-Newtonian (PPN) formalism, could furnish tests for the same theory [U [33, [H, . 

In this paper, we want to face the problem to study /(i?) gravity at a fundamental level. In particular, in the 
framework of the Arnowitt-Deser-Misner {AVM) formalism [42|, we want to investigate the possibility to find out 
cosmological terms as eigenvalues of generalized /(i?)-IIamiltonians in a Sturm-Liouville-likc problem^. This issue is 
particularly relevant from several viewpoints. First of all, our aim is to show that vacuum energy of gravitational 
field is not a particular feature of GR where the cosmological constant has to by added by hand into dynamics. At a 
classical level, it is well known that f{R) gravity, for the Ricci scalar R equal to a constant, exhibit several deSitter 
solutions 01 but a definite discussion, at a fundamental level, considering cosmological terms as eigenvalues of such 
theories is lacking. Besides, the computation of the Casimir energy, the seeking for zero point energy in different 
backgrounds^ give a track to achieve one-loop energy regularization and renormalization for this kind of theories 
(4^ [4]|. On the other hand, these issues can be considered in a multigravity approach to spacetime foam if the A/" 
spacetimes constituting the foam are supposed to evolve, in general, with different curvature laws and ground states 
(cosmological constants) 0]. 

The layout of the paper is the following. In Sec. II, we recall the Hamiltonian formalism in the (AVAi) approach 
of GR. It is developed for generic f{R) gravity in Sec. III. Sec. IV is devoted to find out the cosmological constant as 
the eigenvalue of a generalized f{R) Hamiltonian. We discuss the orthogonal decomposition of the wave functional 
and derive the total one-loop energy density for the transverse-traceless tensor component. In Sec.V, we give an 
example: the transverse-traceless spin 2 operator is calculated for the Schwarzschild metric and the energy density 
contributions to the cosmological constant are calculated in the WKB approximation. This is a realization of the 
above formal cosmological constant calculation. Sec. VI is devoted to the one-loop energy calculation by the zeta 
function regularization method. The explicit value of the renormalized Aq constant, considered as a running constant, 
is achieved. How it can be set to zero is explicitly derived for f{R) = exp{—aR). Summary and conclusions are drawn 
in Sec. VII. In Appendix A, details on zeta function regularization are given. 

II. THE HAMILTONIAN CONSTRAINT OF GENERAL RELATIVITY 

Let us briefly report how to compute the Hamiltonian constraint for GR considering the standard Hilbert-Einstein 
theory / (i?) = R and the Arnowitt-Deser-Misner {ADM ) 3 -I- 1 decomposition ^] . In terms of these variables, the 



^ Following the Occam razor prescriptions: "Entia non sunt multiplicanda praeter necessitatem. " 

^ See Ref. [431, for the application of the Sturm-Liouville problem in the simple case of / (R) = R, even in presence of a massive graviton. 
^ For different f{R), we expect different zero point energies and, obviously, different vacuum states. 
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line element is 

ds^ = (x) dxf'dx'-' = {-N^ + N,N^) dt^ + 2Njdtdx' + g.jdx'dx^ . 

N is the lapse function, while Ni the shift function. In terms of these variables, the gravitational Lagrangian, with 
the boundary terms neglected, can be written as 

£ [N, N,,g,,] = = -^-^ [K^JK'' ~ + {^R - 2A,)] , (1) 

where Kij is the second fundamental form, K = g^^ Kij is the trace, '^R is the three dimensional scalar curvature and 
'^g is the three dimensional determinant of the metric. The conjugate momentum is simply 



By a Legendre transformation, we calculate the Hamiltonian 

H = Jd^x [NH+N,H'] , (3) 

where 



n = {2k) GyHTT^^'^' -^{'R- 2Ae) (4) 

and 

n' = ~2vy\ (5) 

where Aj, is the bare cosmological constant. The equations of motion lead to two classical constraints 

n = o 



W ^0 ' 

representing invariance under time re-parameterization and invariance under diffeomorphism, respectively. Gijki is 
the supermetric defined as 

Gijki = 2^''^*''^''' ^ 9ii9jk - gijgu)- (7) 

When 7i and TV are considered as operators acting on some wave function, we have 

[g,,] = (8) 

and 

H,* [5,,] = 0. (9) 

Eq.® is the Wheeler-De Witt equation ( WDW) (iH] . Eqs.® and Q describe the wave function of the universe ^ [gij]- 
The WDW equation represents invariance under time re-parameterization in an operatorial form. This standard lore 
can be applied to a generic f{R) theory of gravity with the aim to achieve a cosmological term as an eigenvalue of 
the WDW equation. 



III. THE HAMILTONIAN CONSTRAINT FOR A GENERIC f{R) THEORY OF GRAVITY 



Let us consider now the Lagrangian density describing a generic f{R) theory of gravity, namely 

C = ^g{f{R)~2K,), with/VO, 



(10) 
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where / (i?) is an arbitrary smooth function of the scalar curvature and primes denote differentiation with respect 
to the scalar curvature. A cosmological term is added also in this case for the sake of generality. Obviously /" = 
corresponds to GR. The generalized Hamiltonian density for the / (R) theory assumes the form^ 



^ - 2Ae - 3K,jK'^ + K^) + V{V) - ^g^'V^,, - 2f^K, 



where 



V{V)^^[Rf'iR)-f{R)] 



(11) 



(12) 



Henceforth, the superscript 3 indicating the spatial part of the metric will be omitted on the metric itself. When 
/ (i?) = R, V{V) = as it should be. Since 



V'^ = -2v^5'V' [R) =^ V^- 6Vff/' (i?) , 

we have 

With the help of Eq.®, Eq.dll]) becomes 



(13) 



(14) 



n = f (R) 

However 



(2k) GijkiTT^-' tt'^'' 



2k \ 



(3)i?-2A, 



1 

2k, 



Vgf (R) {2K,,K'^) + V{V) + 2g^= {^gf (i?)),,^- - 2p'i K,, 

(15) 



(16) 



then we obtain 



n = f' (R) 



(2K)Gyfe,7r'^'^'=' 



2k V 



(3)i?_2A, 



2k 



2VgK,,K^' (/' (i?) - 1) + V{V) + 2g^^ {Vgf (i?)),, (17) 



and transforming into canonical momenta, one gets 



n = f (R) 



(2«)G.,,z^''-'"7r^'-:^((3)i?-2A 



+2 (2k) 



By imposing the Hamiltonian constraint, we obtain 



f{R) 



{2k)G,,m7T^^7T>^'-^(^^R 



2k 



(/' (R) 1)+^ [v{V) + 2g^^ (Vaf mi., 

(18) 
(19) 



2 (2k) 



(/' (R) - 1) + :^ \viV) + 2g^= iVgf iR))J = -f (R) V9- (20) 



If we assume that /' (R) ^ the previous expression becomes 



Zk 



-(2n) 



2(/'(i?)-l) 



1 



/' [R) 2k/' (i?) 



V{V) + 2g^^{V9r (.R))i 



Now, we integrate over the hypersurface S to obtain 



(2k)G,,,,^^^-^'='-||(3)^ 



(2k) 



TT 2 



Hf{R)-i) 
,/' (R) 



Va—- 

K 

(21) 
(22) 



See also Ref . |46|| for technical details. 
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The term 



(23) 



(24) 



appears to be a three-divergence and therefore wiU not contribute to the computation. The remaining equation 
simphfies into 



Zk 



(2k) 



2(/'(i?)-l) , V{V) 



f (R) 2Kf' (R) 



By a canonical procedure of quantization, we want to obtain the vacuum state of a generic f{R) theory. 



(25) 



IV. THE COSMOLOGICAL CONSTANT AS AN EIGENVALUE FOR THE GENERALIZED / (R) 

HAMILTONIAN 

The standard WDW equation ([5]) can be cast into the form of an eigenvalue equation 

Ae*[5«j] = A(f)*[(7„], 

where 

As = (2k) G,,M7T'^n^' 



2k 



R. 



(26) 



(27) 



If we multiply Eq. ([26|) by 5** [gij] and we functionally integrate over the three spatial metric gij, we get 

V [g,,] [5,,] As* [g,,] ^ jv [g,,] A (x) [g.,,] vJ- [g,,] (28) 
and after integrating over the hypersurface E, one can formally re- write the modified WDW equation as 



1 /P[g„]vE'*[5y]/5.d3a:As*[gy] 1 



Ac 

K 



V JV[g,,]^*[g.,,]^[g.,,] V 
where the explicit expression of A (x) has been used and we have defined the volume of the hypersurface E as 



V = 



(29) 



(30) 



The formal eigenvalue equation ([29|) is a simple manipulation of Eq.®. We can gain more information considering a 
separation of the spatial part of the metric into a background term, gij, and a quantum fluctuation, , 



g,j +h,j. 



Thus Eq.(129l) becomes 



aL°^+a« + aL^) 



(31) 



(32) 



where A^"* represents the i*** order of perturbation in hij. By observing that the kinetic part of As is quadratic in 
the momenta, we only need to expand the three-scalar curvature / d^x^R^^^ up to the quadratic order and we get 



--hAh+-h^'Ahu - i/i'^V,V,/i-- 
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(33) 



where h is the trace of hij and is the three dimensional scalar curvature. By repeating the same procedure for 
the generalized WDW equation Ea.(p5)). we obtain 



1 (* 










V 







') 2n2if'iR)^l){^\J^d^x[G,,ki7r^'7r>^' + f]\^) 1 (* |/s '^'^S^ | *) _ A, 



V /' (i?) 



V 



From Ea. ([M|) . we can define a ^^modified^^ A^' operator which includes /' (_R). Thus, we obtain 



A 



(2) 



^ /' (^) 



(M-lvP) 



where 



with 



A 



(2) 



(2K)/i(i?) G,yM^*-''7r^'' 



2k 



h{R)^l 



2[r(i?)-i] 



(34) 



(35) 



(36) 



(37) 



and where ^/j^™ jg the linearized scalar curvature whose expression is shown in square brackets of Eg. ((33|) . Note that 
when / (i?) = R, consistently it is h (R) — 1. From Eg. ([55]) . we redefine Af, 



Al = Ae + 



1 



r d'^x^ 



2V 



2VJ^ f'iR) 



(38) 



where we have explicitly used the definition of V(V). In order to make explicit calculations, we need an orthogonal 
decomposition for both 7r,jj and hij to disentangle gauge modes from physical deformations. We define the inner 
product 



{h,k):= / ^G'^'''h,,{x)kki{x)d^x, (39) 

by means of the inverse WDW metric Gijki , to have a metric on the space of deformations, i.e. a quadratic form on 
the tangent space at hij, with 



Qijkl ^ ^(g^kgjl ^ gUgjk „ 2 ^ ) . 

The inverse metric is defined on cotangent space and it assumes the form 



(40) 



(P, q) ■■= / ^/gG^JklP"' {x) q^^ {x) d^x, 



(41) 



SO that 



(42) 



Note that in this scheme the "inverse metric" is actually the WDW metric defined on phase space. The desired 
decomposition on the tangent space of 3-metric deformations (i^, IH, [sOl is: 



1 



hg,, + (LO 



(43) 
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where the operator L maps into symmetric tracefree tensors 



= V,e, + - ^9^, (V • ■ 



Thus the inner product between three-geometries becomes 



(44) 



V5 



With the orthogonal decomposition in hand we can define the trial wave functional as 



where 



* [h, (1?)] = (1?)] VP 4 (1?) VP [/i*™- (^)] , 



vl/[^.(l?)] =exp{-i(/ii^-i/i)^J 
VP (1?)] = exp {-i {hK'^hf-"} 



(45) 



(46) 



(47) 



The symbol "_L" denotes the transverse-traceless tensor (TT) (spin 2) of the perturbation, while the symbol "||" 
denotes the longitudinal part (spin 1) of the perturbation. Finally, the symbol "irace" denotes the scalar part of the 
perturbation. A/" is a normalization factor, (•, •)^ denotes space integration and is the inverse "propagator^\ We 
will fix our attention to the TT tensor sector of tne perturbation representing the graviton. Therefore, representation 
(|^ reduces to 



■H [hij (^)] = A/'exp <^ -- (hK-^h) 



(48) 



Actually there is no reason to neglect longitudinal and trace perturbations. However, following the analysis of 
Refs.[49i. i51., ,52j on the perturbation decomposition, we can discover that the relevant components can be restricted 
to the TT modes and to the trace modes. Moreover, for certain backgrounds, TT tensors can be a source of instability 
as shown in Refs.[5l|, [H, H^. Even the trace part can be regarded as a source of instability. Indeed this is usually 
termed conformal instability. The appearance of an instability on the TT modes is known as non-conformal instability. 
This means that does not exist a gauge choice that can eliminate negative modes. To proceed with Eq. (|32|) . we need 
to know the action of some basic operators on vp [hij]. The action of the operator hij on |vp) = vp [ft,,-^] is realized by 

M 



%(a;)|*) = /iy(^)^'[/i„-]. 



(49) 



The action of the operator ntj on |vp), in general, is 



TT,, (a;)!*) 



while the inner product is defined by the functional integration: 



We demand that 



V 



V 



/I?[5y]vp* [%]vp[/i,,] 



(50) 



(51) 



(52) 
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be stationary against arbitrary variations of [hij]. Note that Eq. (|52|) can be considered as the variational analog 
of a Sturm-Liouville problem with the cosmological constant regarded as the associated eigenvalue. Therefore the 



solution of Eg. ([^5]) corresponds to the minimum of Ea.((52|). The form of ( ^E" 
of the wave functional ([48l) and with the help of 



A. 



and 



{^\h,ji^)hki (^)l^) 



= 



^ ) can be computed with the help 



(53) 



Extracting the TT tensor contribution, we get 



The propagator {x,x)^^i^i can be represented as 



(^,^) 



iakl 



E 



2A(r) 



(54) 



(55) 



(56) 



where h^iJ^ {~x) are the eigenfunctions of A2. t denotes a complete set of indices and A (r) are a set of variational 
parameters to be determined by the minimization of Eq.(|55p. The expectation value of is easily obtained by 
inserting the form of the propagator into Ea. (|55p 



A^AO 1 



EE 



(2K)/i(i?) A, (r) 



(2k) a, (t) 



(57) 



By minimizing with respect to the variational function Ai (t), we obtain the total one loop energy density for TT 
tensors 



A^ (A,) = -n^/h{R)-J2 + 



(58) 



where A^ is expressed by the Eq.(l38|). The above expression makes sense only for ujf (r) > 0. It is the main formal 
result of this paper. It is true for generic /(i?) functions since h{R) explicitly appears in it. 



V. THE TRANSVERSE TRACELESS (TT) SPIN 2 OPERATOR FOR THE SCHWARZSCHILD 

METRIC AND THE WKB APPROXIMATION 

The above considerations can be specified choosing a given metric. For example, the quantity A^ can be calculated 
for a Schwarzschild metric in the WKB approximation. Apparently, there is no a strong motivation to consider a 
Schwarzschild metric as a probe for a cosmological problem. Nevertheless, every quantum field induces a "cosmological 
term" by means of vacuum expectation values and the variational approach we have considered is particularly easy 
to use for a spherically symmetric metric. The Schwarzschild metric is the simplest sourceless solution of the Einstein 
field equations which can be used to compute a cosmological constant spectrum. Of course, also Minkowski space 
can be put in the form of a spherically symmetric metric, but in that case there is no gravity at all. The other 
solution need a source which is not considered in the present paper. In this sense, the computation is a real vacuum 
contribution to the cosmological term. The spin-two operator for the Schwarzschild metric is defined by 

(A^/i^^)^ := - [t^Tk'^t + 2 {Rh'^y^ , (59) 
where the transverse-traceless (TT) tensor for the quantum fluctuation is obtained by the following decomposition 

H^H- \5lh + \5lh = {h^)\ + ^Ih. (60) 
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This implies that {h^)^ (5* = 0. The transversahty condition is apphed on {h^Y^ and becomes Vj {h^Y^ ~ 0. Thus 

- (Arh^-Y = -As {h-ll + 5 fl - ^) - (61) 



where A5 is the scalar curved Laplacian, whose form is 

2MG\ (P 



A. = 1 - 



and R°: is the mixed Ricci tensor whose components are: 



2r - ZMG 



dr 



Rf = <- 



2MG MG MGj 



(62) 



(63) 



This implies that the scalar curvature is traceless. We are therefore led to study the following eigenvalue equation 

{A,h^^Y=^'^) (64) 

where to^ is the eigenvalue of the corresponding equation. In doing so, we follow Regge and Wheeler in analyzing 
the equation as modes of definite frequency, angular momentum and parity 55|. In particular, our choice for the 
three-dimensional gravitational perturbation is represented by its even-parity form 



with 



(/i™^")!. (r, 0) = diag [H (r) ,K{r),L (r)] Yi^ {d, , 



H{T)^h\ (r)-i/i(r) 
K{r)^hl (r)-k(r-) 
L{r)^hl{r)-lh{r) 



(65) 



(66) 



From the transversahty condition, we obtain h\ (r) = h\ (r) . Then K (r) ^ L (r) . For a generic value of the angular 
momentum L, representation ([65]) joined to Eq. ([6T|) lead to the following system of PDE's 



(-A5 + ^ 

f-Ac + St 



(67) 



Defining the "reduced" fields 



H{r) = 



Mr) 



K{r) 



h{r) 



r r 
and passing to the proper geodesic distance from the throat of the bridge 

dx — ± 



1 



2MG 



(68) 



(69) 



the system ([67)1 becomes 



with 



-i^ + Vi (r)J h [x) - u^l J, {x) 
~£ + V2 (r)l h ix) = co^h ix) 



Vi (r) = ^ + C/i (r) + m 
V2 (r) = ^ + U2 (r) + TO 



(70) 



(71) 
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where we have defined r = r (x) and 



(72) 

f/2(r)=[^(l-2MG) + 3MG] 



Note that 



Ui (r) > when r > 
Ui (r) < when 2MG <r < ^ ^.^3^ 

U2 (r) > Vr [2MG, +00) 

In order to use the WKB approximation, we define two r-dependent radial wave numbers fci (r, Z, wi ^i) and 

kf (r, I, uJi^ni) = i^j.ni - - ml (r) 

(74) 

where we have defined two r-dependent effective masses (r) and (r). The WKB approximation we will use to 
evaluate Eq. (l58p is equivalent to the scattering phase shift method and to the entropy computation in the brick wall 
model. We begin by counting the number of modes with frequency less than cui, i — 1,2. This is given approximately 
by 



9i^^) = J ly^{l,L^^){2l + l), (75) 

where j/^ {I, uji), i = 1, 2 is the number of nodes in the mode with (Z, uji), such that (r = r (x)) 

1 /■+°° I 

Vi{l-,^i) = 7r dx\kf{r,l,iOi). (76) 

Here it is understood that the integration with respect to x and I is taken over those values which satisfy fc| {r,l,u!i) > 0, 
i = 1,2. With the help of Eqs. (|75| 76p . we obtain the one loop total energy for TT tensors which is 



1 ^ 

i— 1 



2 ^ + cx> 

dx 



+ 00 







dg (wi) 

LU., ; dhJi 



(77) 



(78) 



By extracting the energy density contributing to the cosmological constant, we get 

^ f r+oo . /•+00 . 

K = K,i + K,2 = + p2 = -v/HR)y^ _ ^2 ^^^^^^ _^ ^A^u^l- (r)du;2 

where we have included an additional 47r coming from the angular integration. 

VI. ONE-LOOP ENERGY REGULARIZATION AND RENORMALIZATION 

In this section, we will use the zeta function regularization method to compute the energy densities p\ and p2- 
Note that this procedure is completely equivalent to the subtraction procedure of the Casimir energy computation 
where the zero point energy (ZPE) in different backgrounds with the same asymptotic properties is involved. To this 
purpose, we introduce the additional mass parameter /i in order to restore the correct dimension for the regularized 
quantities. Such an arbitrary mass scale emerges unavoidably in any regularization scheme. Then we have 

K /■+°° i.P- 



(e) = -^h{R) \ dco, ^ (79) 

167r2 Jo (^2„^2(^))^-2 
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where 



(80) 



The integration has to be meant in the range where ujf — mf (r) > 0^. One gets 



ffl^ ( 



+ 1^ ( 2 / \ 



2 In 2 - i 
2 



(81) 



i = 1,2. In order to renormahze the divergent ZPE, we write 

= SttG [pi (e) + P2 (e) + Pi (m) + P2 (m)] , (82) 

where we have separated the divergent part from the finite part. For practical purposes, it is useful to divide with 
the factor y/h (R). To this aim, we define 

a: 



Rf {R) - f (R) 



and we extract the divergent part of A, in the limit e — > 0, by setting 

G 



1 



A' 



div 



SttG [pi (e) + p2 (e)] 



327re 



[mf (r) + 7712 ("r)] ■ 



(83) 



(84) 



Thus, the renormalization is performed via the absorption of the divergent part into the re-definition of the bare 
classical cosmological constant Ac, that is 



Ac Ao + VHR)A''™. 
The remaining finite value for the cosmological constant reads^ 



(85) 



KM 

SttG 



Pi M + P2 (p) 



2567r^ 



777^ (r) 



In 



|777f (r)| 



2 In 2 



where 



-7772 ('') 



In 



777^ (r) 



2 In 2 - - 
2 



K (m) 



A / N 1 /■ ,3 ^Rf'{R)-f{R) 



f'iR) 



(86) 



(87) 



is the modified cosmological constant. The quantity in Eq. ([5S)) depends on the arbitrary mass scale /i. It is app ropriate 
to use the renormalization group equation to eliminate such a dependence. To this aim, we impose that [56[ 

Solving it, we find that the renormalized constant Aq should be treated as a running one in the sense that it varies, 
provided that the scale /z is changing 



G 



A^ (^, r) = A'o {^lo, r) + ^ [mf (r) + (r)] In 

Idtt pq 



(89) 



^ Details of the calculation can be found in the Appendix A. 

^ Since rrij (r) can change in sign, when we integrate over uii we can use either /-|- or /_ . This leads to the appearance of the absolute 
value. 
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Substituting Eq.(189l) into Eq.jMl) we find 



SttG 



2567r2 



rrii (r) 



In 



|™i (^) I 
M5 



- 2 In 2 



- (r) 



In 



"ij jr) 
f4 



- 2 In 2 - 



(90) 



It is worth remarking that while rn^ (r) is constant in sign, mf (r) is not. Indeed, for the critical value r — 5MG/2, 
ml (f) = rrig and in the range {2MG, 5MG/2) for some values of rrig, m\ (f) can be negative. It is interesting therefore 
concentrate in this range. To further proceed, we observe that m\ (r) and (r) can be recast into a more suggestive 
and useful form, namely 



ml (r) = Ui (r) = (r, M) - ml (r, M) 
ml (r) = U2 (r) = (r, M) + ml (r, M) 



(91) 



where m\ {r, M) when r ^ 00 or r — > 2MG and 7712 {r, M) — iMG/r^. Nevertheless, in the above mentioned 
range m\ (r, M) is negligible when compared with ml (r, M). So, in a first approximation we can write 



nl (r)~-m2 {ra,M) = -ml (M) 
ml (r) ~ ml {ro,M) = TOq (M) 



(92) 



where we have defined a parameter ro > 2MG and m§ (M) = SMG/rj^. The main reason for introducing a new 
parameter resides in the fluctuation of the horizon that forbids any kind of approach. Of course the quantum 
fluctuation must obey the uncertainty relations. Thus Eq. ([M)) becomes 



Af,(Mo,r)_ mtiM) 



In 



SttG 1287r2 
Now, we compute the maximum of Ag, by setting 

_ ml (M) 



ml (M) 



Thus An becomes 



Ag {fJ.o,x) 



In (x) 



(93) 



(94) 



(95) 



As a function of x, Aq (/xq, x) vanishes for a; = and x — exp (— ^) and when x [O, exp (—5)] , Ag {fiQ, x) > 0. It has 
a maximum for 



1 



ml (M) = ^ 



and its value is 

Ao(/io,.x) = ^ 
Isolating Ag (/ig, x), we get 



VHR) 



Ag {^io,x) 



1 



dPx^ 



Rf (R) - / (R) 
f'iR) 



27re2 



Ao(Mo,x) = ^AMg|-^ 



d?x^ 



Rf (R) - f (R) 



Note that Ag (/ig, x) can be set to zero when 

■ g4 



27re2 



— / (fx^ 
2Vj^ 



RfiR)-fiR) 
f'{R) 



(96) 



(97) 



(98) 



(99) 



13 



Let us see what happens when / (i?) = exp {—aR). This choice is simply suggested by the regularity of the function 
at every scale and by the fact that any power of R, considered as a correction to GR, is included. In this 
case, Eq.(j99| becomes 



/3a exp {-aR) + 2 G/ij^ _ 1 f ^ 



y Q!exp(— ai?) Tre^ 
For Schwarzschild, it is i? = 0, then 



^7 / d^Xy^{l + aR). (100) 



V(3a + 2)a- ^^^^^ 



By setting a — G, we have the relation 



^"-^V(3GT2)G- (10^) 
Remark Note that in any case, the maximum of A corresponds to the minimum of the energy density. 



VII. SUMMARY AND CONCLUSIONS 



Despite of the successes of General Relativity, such a theory can only be considered as a step toward a much more 
complete and comprehensive structure due to a large number of weaknesses. Among them, the issue to find out the 
fundamental gravitational vacuum state is one of the main problem to achieve a definite Quantum Gravity theory 
which, till now is lacking. However, several semiclassical approaches have been proposed and, from several points 
of view, it is clear that the former Hilbert-Einstein scheme has to be enlarged. The f{R) theories of gravity are a 
minimal but well founded extension of GR where the form of the function f {R) is not supposed "a priori" but is 
reconstructed by the observed dynamics at galactic and cosmological scales [3l|, [M] . Also if they seems a viable scheme 
from cosmology and astrophysics viewpoints, their theoretical foundation has to be sought at a fundamental level. 
In particular, one has to face the possibility to encompass the f{R) gravity in the general framework of Quantum 
Field Theory on curved spacetime. In this paper, we have dealt with the problem to find out vacuum states for f{R) 
gravity via the (3 + 1) ADM formalism. Analogously to GR, we have constructed the Hamiltonian constraint of a 
generic f{R) theory and then achieved a canonical quantization giving the /(i?)-WDW equation. In this context, the 
cosmological constant (vacuum state) emerges as a WDW eigenvalue. The related wave functional can be split by an 
orthogonal decomposition and then, constructing the transverse-traceless propagator, it is possible to obtain, after a 
variational minimization, the total one-loop energy density for the TT tensors. Such a quantity explicitly depends on 
the form of f{R)- As an application, we derive the energy density contributions to the cosmological constant for a 
TT spin 2 operator in the Schwarzschild metric and in the WKB approximation. The one-loop energy regularization 
and renormalization are achieved by the zeta function regularization method. The resulting renormalized Ag is a 
running constant which can be set to zero depending on the value of an arbitrary mass scale parameter /i. As explicit 
calculation, we find out the value of such a parameter for a theory of the form f{R) = exp(— ai?) in the Schwarzschild 
metric. This case can be used for several applications at cosmological and astrophysical scales. In 
particular, truncated version s of such an exponential function, power law f{R), have been used for 
galactic dynamics [32l. l34l |35^. In those cases, a corrected Newtonian potential, derived from the f{R) 
Schwarzschild solution, has been used to fit, with great accuracy, data from low surface brightness 
galaxies without using dark matter haloes. This approach has allowed to fix a suitable mass scale 
comparable with the core size of galactic systems. Such a mass can be directly related to the above 
parameter a depending on the core radius Tc (see also [35]). 

In summary, the application of Quantum Field Theory methods to f{R) gravity seems a viable scheme and gives 
positive results toward the issue to select vacuum states (eigenvalues) which can be interpreted as the cosmological 
constant. However, further studies are needed in order to generalize such results to other metrics and other Extended 
Theories of Gravity. 



APPENDIX A: ZETA FUNCTION REGULARIZATION 
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In this appendix, we report details on computation leading to expression (|79[) . We begin with the following integral 

1+ = f^'^" lo^ duj- 



with (r) > 0. 



1. J+ computation 



If we define t — LLij ^rn? (r), the integral in Eg. ljAip becomes 



/' + OC 



p(e)= ii'^m^-'^ (r) 



(t2 + ir 



1 2e 4-2S , , r(i)r(£-2) 

-u m (r) ; 



(Al) 



,2 \ £ 



where we have used the following identities involving the beta function 



Jo 

related to the gamma function by means of 



+ 00 



dt 



B{x,y) 



V{x)T{y) 



Rex > 0, Rey > 



r (.T + y) 

Taking into account the following relations for the F-function 

Vie -2)- 



and the expansion for small e 



r (1 + e) = 1 - 7e + O (e^) 

r(£+^)= r(i)-er(i)(7 + 21n2) + 0(e2), 
= 1 + elna; + O (e^) 



where 7 is the Euler's constant, we find 



p{e) = - 



TO* (r) 
16 



£ [r) 



2 In 2 - 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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2. 7 computation 



If we define t = uj/ yjw? (r), the integral /_ in Eq. ljAip becomes 



p{e) = ^l''m^-'' {r) / dt 



+ CXD 



(t» 



where we have used the following identity involving the beta function 

is (i - - ^) = /+°° dtt^-^ {tp - ly-' 

p > 0, Rei^ > 0,Re^ < p — pRev 

and the reflection formula 

r (z) r (1 - z) = -zT (-z) r (z) 

From the first of Eqs. (jA5p and from the expansion for small e 

r(^-e) =r0) (l-e(-7-21n2 + 2)) + 0(£2 



(A8) 



(A9) 



(AlO) 



1 + elnx + O (e^) , 



we find 



Pie)-- 



(r) 
16 



+ 1 2 ^ ^ 



2 In 2 - - 
2 



(All) 



(A12) 
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